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Abstract. We report on high-resolution spectroscopy of ultracold fermionic 
^^Na'^^K Feshbach molecules, and identify a two-photon pathway to the rovibrational 
singlet ground state via a resonantly mixed intermediate state. 

Photoassociation in a ^^Na-^°K atomic mixture and one-photon spectroscopy on 
^^Na'^^K Feshbach molecules reveal about 20 vibrational levels of the electronically 
excited c^E+ state. Two of these levels are found to be strongly perturbed by 
nearby B^II states via spin-orbit coupling, resulting in additional lines of dominant 
singlet character in the perturbed complex B^n|u=4)~c^E+|?;=25), or of resonantly 
mixed character in B^n|u=12)~c^E+|u=35). The dominantly singlet level is used 
to locate the absolute rovibrational singlet ground state X^E+|w=0, J=0) via Autler- 
Townes spectroscopy. We demonstrate coherent two-photon coupling via dark state 
spectroscopy between the predominantly triplet Feshbach molecular state and the 
singlet ground state. Its binding energy is measured to be 5212.0447(1) cm“^, a 
thousand-fold improvement in accuracy compared to previous determinations. In their 
absolute singlet ground state, molecules are chemically stable under binary 

collisions and possess a large electric dipole moment of 2.72 Debye. Our work thus 
paves the way towards the creation of strongly dipolar Fermi gases of NaK molecules. 
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1. Introduction 


Ever since laser cooling and evaporative cooling gave full control over the motional 
and internal degrees of atoms, there has been a strong effort to extend such control 
over the richer internal structure of molecules. With large samples of molecules, all 
occupying the same internal quantum state, controlled switching of chemical reactions 
could be studied at the quantum level. Their internal degrees of freedom make 
molecules suitable as potential carriers for quantum information, as a quantum resource 
for precision measurements and as the building block for novel quantum many-body 
systems. A particular appeal lies in the creation of new forms of quantum dipolar matter, 
such as topological superfluids or quantum crystals. The enormous progress towards 
ultracold molecules over the last decade has been summarized in review papers HI. 
to which one may add the recent successes of magneto-optical trapping [^[^, a novel 
type of Sisyphus cooling [^, as well as evaporative cooling of molecules [^, among 
others. An alternative to direct cooling of molecules, which has so far been limited to 
temperatures of several millikelvin, is the creation of weakly bound Feshbach molecules 
at nanokelvin temperatures from a gas of ultracold atoms, which has been achieved for 
homonuclear |8 14 as well as for heteronuclear molecules 15-22 . Feshbach molecules 
created out of fermionic atoms were remarkably long-lived |^, allowing the observation 
of Bose-Einstein condensation 23 28 . However, the dipole moment of such long-range 
Feshbach molecules is still vanishingly small. In pioneering works at JILA and in 
Innsbruck, ultracold molecules of fermionic “^^K^^Rb and bosonic non-polar ^^^Cs 2 in 
their absolute rovibrational ground state were created from weakly bound Feshbach 
molecules via a coherent two-photon transfer 29 31 . ^°K®^Rb molecules were found to 
undergo the reaction 2KRb —)■ K 2 -|- Rb 2 32 , precluding further cooling into the deeply 


degenerate regime. However, conhnement to two- and three-dimensional optical lattice 


potentials allowed to stabilize molecular gases of KRb against reactive collisions 33,34 


Recently, ultracold gases of bosonic RbCs have been created 35,36 , a molecule that is 


chemically stable under binary collisions 37 


In the present work we focus on the fermionic molecule ^^Na‘^°K, also chemically 


stable under two-body collisions 37 and known to possess a large electric dipole moment 
of 2.72 Debye 38 , about hve times larger than that of KRb. It is therefore an ideal 


candidate for the formation of a strongly interacting, stable dipolar Fermi gas with 
interaction energies on the order of tens of percent of the Fermi energy. Following our 
creation of ultracold Feshbach molecules 20 , we here perform one- and two-photon 


spectroscopy on ^^Na^°K and identify a two-photon pathway from the predominantly 
triplet Feshbach molecular state to the absolute rovibrational singlet ground state. A 
priori, the existence of such a pathway is not obvious, as spin-orbit coupling, mixing 
singlet and triplet states, is weak for this light molecule. 

Historically, NaK was a textbook example of a diatomic molecule and one of the 
best studied diatomic molecules before the advent of photoassociation and coherent two- 
photon spectroscopy of cold and ultracold KRb. The earliest absorption spectra 39 
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predate modern quantum mechanics, but already within a decade the singlet ground and 
first excited Born-Oppenheimer potentials had been experimentally characterized 40 


In particular, levels of the B^II state were already known to an accuracy of 1 cm“^. 
In the era of laser spectroscopy, this accuracy improved to 0.001 cm“^, resulting 
in highly precise excited singlet potentials, for example of the B^fl state relevant 
here 41,42 , and of the singlet and triplet ground states, and a^E’*' 43,44 


In addition, singlet-triplet perturbations were observed in electronically excited states 
that led to fluorescence down to the triplet ground state, first D^fl —)• a^E+ 45 -48 , 
later B^fl —a^E"*" 41,49-53 . This already demonstrates that two-photon coupling 


is possible between the triplet and singlet ground states, and indeed, one of the 
perturbations studied here (albeit in the previously unstudied isotopologue ^^Na"^°K) was 
first observed in laser-induced fluorescence in [^. The excited triplet states that caused 
such perturbations were characterized via deperturbation techniques (^, leading to 


deperturbed potential energy curves for the b^fl state 55 and the c^E+ state |4ll|56||57|. 


These potentials and calculated energy levels are used in the present study as a guide 
for our spectroscopic search for perturbations in ^^Na^°K, which has not been studied 
before due to ^°K’s low natural abundance of 0.01%. Recently, two-photon spectroscopy 
in the “traditional direction” from the absolute singlet to low-lying triplet ground states 


of NaK was performed on ^^Na^®K in a molecular beam 58 


We here demonstrate that even the highest vibrationally excited triplet ground 
state near a Feshbach resonance can be coherently coupled to the absolute rovibrational 
singlet ground state with significant coupling strength, enabling the formation of ground 
state dipolar molecules from ultracold Feshbach molecules. The potential energy curves 
are shown in Fig. The explored pathway connects the predominantly triplet Feshbach 
molecular state to a perturbed level of mixed B^n~c^E+ character via the up-leg laser 
of frequency Ui. From there, the down-leg laser of frequency U 2 promotes the coupling 
to the absolute singlet ground state. 

The starting point of the experiment is a near degenerate Bose-Fermi mixture of 
23Na-^°K, as described in our previous work (Ml. Since in heatpipe or beam experiments 


only the deep singlet molecular levels are populated, the triplet states c^E+ had only 
been observed where they were strongly perturbed by a nearby B^fl state. There are 
thus uncertainties in the exact location of c^E"*" energy levels of more than 1 cm“^. The 
first task is therefore to perform photoassociation spectroscopy on the atomic ultracold 
Na-K mixture, as described in section This allows the direct observation of c^E"*" 
levels, and thanks to the ultracold temperatures, only low lying rotational states are 
accessed. Equipped with the exact locations of the c^E+ levels, and highly accurate 
predictions on the locations of singlet B^fl levels from mass-scaled Dunham coefficients 
for ^^Na^^K 41,42 , we identify two accidental degeneracies between c^E"*" and B^fl in 
^^Na^°K, namely B^n|u=4)~c^E+|u=25) and B^n|u=12)~c^E+|u=35), where v denotes 
the vibrational quantum number. In section we therefore perform high-resolution 
spectroscopy on these near-degenerate levels starting with Feshbach molecules. We 
resolve detailed fine- and hyperfine structure that is well reproduced by our theoretical 
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Figure 1. Ground and excited state molecular potentials of NaK relevant to the two- 
photon STIRAP transfer to the absolute rovibrational ground state |G) in the X^I]+ 
potential. The initial state of weakly bound Feshbach molecules |F) is associated 
with the highest vibrational level |u=19) of the a^E+ potential. Excited states |E) 
explored in this work are dominantly associated with the c^E+ and B^II potentials, 
and also influenced by states of the b^II potential. As an example, the vibrational 
wavefunctions of the coupled complex B^n|u=4)~c^E+|w=25) are shown. The states 
of the three-level system |F), |E) and |G) are coupled by up-leg and down-leg lasers 
with frequencies wi, and Rabi couplings fli, ^ 2 - 


calculations. The manifold B^n|u=4)~c^S’''|u=25) is sufficiently perturbed so that the 
singlet-dominant feature can be directly addressed from triplet Feshbach molecules. 
For B^n|u=12)~c^S+|u=35), the observation of two pairs of lines with nearly identical 
hyperhne structure signals exceptionally strong mixing, originating from an essentially 
perfect degeneracy and sufficiently strong spin-orbit coupling. Finally in section 
we perform two-photon spectroscopy, first to locate the absolute rovibrational singlet 
ground state using Autler-Townes spectroscopy, then to perform coherent dark state 
spectroscopy for a high-precision measurement of the binding energy of the rovibrational 
ground state. The strong singlet-triplet mixing in the intermediate state will enable 
stimulated rapid adiabatic passage (STIRAP) from Feshbach molecules to absolute 
singlet ground state molecules. Our work thus opens the door to the creation of a 
strongly interacting dipolar Fermi gas of ^^Na^^K molecules. 
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2. Photoassociation spectroscopy of ^^Na - '^'^K 


401 


Photoassociation (PA) spectroscopy is commonly applied to laser-cooled atomic gases 
trapped in magneto-optical traps (MOTs). However, atomic densities in MOTs are 
limited to about 10^^ cm“^, and free-bound transitions to excited molecular levels can 
typically only reach the long-range part of the excited state potentials. PA in a MOT has 
additional complications for heteronuclear molecules: The excited state is a strong Cq 
potential, which has a shorter range than the potential in homonuclear molecules; and 
overlap of the two MOTs results in severe light-assisted collisions, strongly reducing the 
effective pair density. Instead of using a MOT, we therefore opted for creating a dense 
Bose-Fermi mixture of in an optical trap, cooled close to quantum degeneracy. 

This allows us to work with atomic densities of up to lO^^cm”^. The PA laser can 
illuminate the atomic sample for 10 seconds at high laser power, revealing even deeply 
bound lines of the excited state potentials. 

To summarize the sequence, as outlined in references 20,59 , laser cooled ^^Na 
and atoms are first captured in a magneto-optical trap and optically pumped to 
the |/Na,^/Na) = |2,2) and = |9/2,9/2) states before being transferred into 

an optically plugged magnetic quadrupole trap. In the magnetic trap, the mixture is 
cooled by forced radiofrequency (rf) evaporation of ^^Na on the |2, 2) —)■ |1,1) transition. 

atoms are sympathetically cooled via collisions with ^^Na atoms. At a temperature 
of about 2 /iK, the mixture is loaded into a crossed optical dipole trap operating at 
1064 nm and each species is transferred into the lowest hyperfine states |1,1) for Na 
and |9/2, —9/2) for K via Landau-Zener rf sweeps. For the measurements presented 
in this work, further evaporation in the dipole trap lowers the temperature to about 
1 pK, creating an equal mixture of about 10^ atoms of both species close to quantum 
degeneracy. 

PA spectroscopy is performed on the atomic ^^Na-^°K mixture to explore the 
electronically excited potentials of ^^Na^°K (see Fig.[^. To this end, the near degenerate 
mixture of Na and K is illuminated with an intense laser beam (laser 1), generated by 
a tunable continuous wave titanium sapphire laser. The typical peak intensity at the 
position of the atoms is 7 kW/cm^. The polarization of the light is fixed to be diagonal 
with respect to the vertical magnetic bias field to address both tt and a~^^~ transitions 
to excited states of ^^Na^°K. In order to efficiently cover a large spectral range during 
a single experimental run, the frequency of the PA laser is swept while illuminating 
the ultracold mixture. In general, as we address deeper lying vibrational levels of the 
electronically excited states, the sweep time is extended to overcome the reduction in 
Franck-Condon factors. Typical sweep rates range between 1 GHz/sec to 10 GHz/sec. 
The center positions of the sweeps are recorded using a commercial wavemeter with a 
spectral resolution of 10 MHz and an absolute accuracy of about 300 MHz. 

Photoassociation leads to a simultaneous loss of ^^Na and "‘'’K atoms from the dipole 
trap. After exposure, photoassociation resonances are located by counting the numbers 
of remaining Na and K atoms on absorption images of both species. The cycle time for 
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Figure 2. Photoassociation spectra accessing excited molecular states of ^^Na^°K. 
The observed lines (upper panel) are compared to the prediction of vibrational levels 
(lower panel) of B^II (blue), c^E+ (red), b^IIo (green), with lighter colors denoting 
predictions that lie outside the validity of known Dunham coefficients. Dashed vertical 
lines indicate the position of resonances that are assigned to the c^E+ potential. The 
dot-dashed lines on the right show the Di and D 2 transitions at 12985.2 cm“^ 
and 13042.9 cm“^, respectively. 


a single set of absorption images is about 30 seconds. 

We investigated electronically excited molecular states in the energy range from 
11900 cm“^ up to 13000 cm“^, below the atomic Di transition of Figure 
shows a map of the recorded photoassociation spectra. In order to assign electronic 
states and vibrational quantum numbers to observed PA transitions, the measured 
transition frequencies are compared to predictions based on Dunham coefficients for 
23Na39K [^[^1^ that we mass-scaled for ^^Na'^^K. 

We remind here that previous information on the location of triplet energy levels 
was available only through deperturbation of strongly perturbed levels mixing with 
singlet states. Errors beyond 1 cm“^ are typical, as we see a posteriori from our 
study. Using an ultracold mixture of unbound atoms, we instead directly access triplet 
states. Pairs of ^^Na and atoms in the |1,1) and |9/2, —9/2) states combine into 
a predominantly “stretched” total spin state \F,mp) = \1 /2, —7/2). Unbound atom 
pairs of ^^Na and thus have a dominant (87.5 %) triplet character. Accordingly, 
most of the observed photoassociation resonances can be assigned to the c^S’*' potential. 
This becomes obvious by comparison to predicted positions of vibrational states of the 
relevant excited state potentials b^II, c^S’*', and B^II (see Fig. §. The b^n state is 
expected to be less visible, as the transition dipole moment between a^S’*' and b^II 
vanishes at short-range 60 . From the measured c^S+ energy levels, it will be possible 


in future work to generate a highly accurate potential energy curve for this electronic 
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state. 

The photoassociation spectra allow us to identify two candidate states near 360 THz 
and 372.5 THz, in which vibrational levels of the c^S’*' and B^H potentials are close 
enough to display signihcant singlet-triplet mixing (see Fig. [^. The pairs in question 
are B^n|n=4)~c^E+|n=25) and B^n|n=12)~c^S+|n=35). In both cases, accidentally 
a third electronic state, b^n^, intervenes, where the sub-states = 0,1,2 are split by 
15 cm“^ 55 due to spin-orbit interaction. Here, fl denotes the total orbital and spin 


angular momentum projection along the internuclear axis. 

For strong two-photon coupling, we require an intermediate state that strongly 
connects the initial Feshbach molecular state, dominantly associated with the triplet 
state a^S+|n=19, A^=0, J=l), to the absolute rovibrational singlet ground state 
X^S'''|n=0, J=0). Here, N denotes the quantum number associated with the rotational 
angular momentum, and J the total angular momentum neglecting nuclear spins. There 
are two criteria that need to be fulhlled: First, the intermediate state must feature strong 
singlet-triplet mixing to mediate coupling between the dominantly triplet Feshbach 
molecular state and the singlet rovibrational ground state. Second, the intermediate 
state should have a large transition dipole moment for coupling with both the Feshbach 
and the absolute ground state. Since the electronic part of the transition dipole moment 
is fairly constant for both the a^S+—>■ c^S+ (;^10D) and the B^H—(;^7D) 
transitions 60 , the relevant hgure of merit for the second criterion is the Franck-Condon 
(FC) overlap between the states. Fig. shows the calculated Franck-Condon factors for 
the transitions from a^S+|n=19) to c^S+ and B^H to X^S’''|n=0). The FC factors 
for both candidate states are signihcant. Note that the typical FC factor for the up-leg 
transition is two orders of magnitude smaller than for the down-leg transition, indicating 
that it will be more difficult to achieve substantial coupling between the Feshbach state 
and the excited state than between the excited state and the absolute ground state. 

Regarding the criterion of strong singlet-triplet mixing, spin-orbit matrix elements 
^Bc have been experimentally obtained for B^n~c^S’'' with ^^Na^^K 52,57. These 


works consistently found for B^n|n=4)~c^E’''|n=25) a value of ^bc = 0.16 cm“^ 
and for B^n|n=12)~c^S+|n=35) ^bc = 0.58 cm“^. The spin-orbit matrix elements 
are approximately given by the electronic part of the spin-orbit coupling, which is 
fairly constant as a function of internuclear distance [^, times the Franck-Condon 
overlap between the perturbing states. We checked that this overlap does not change 
signihcantly between ^^Na^^K and ^^Na^°K. Therefore, these values of .^bc give an 
indication how close in energy the singlet and triplet states have to be to undergo 
signihcant mixing. In the next section, we investigate the strength of spin-orbit coupling 
for the candidate states above. 
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Figure 3. Franck-Condon factors between (a) the Feshbach level a^E+|t!=19) and 
excited vibrational levels of the c^E+ potential, denoted by the transition frequency 
wi, and (b) excited vibrational levels of the B^II potential and the absolute ground 
state X^E+|r;=0), denoted by the transition frequency uj 2 - The Franck-Condon factor 
is defined as the squared absolute value of the wavefunction overlap. The vibrational 
wavefunctions are calculated for ^'^Na'^°K using the molecular potentials known for 
23Na39K. The dashed and dot-dashed lines correspond to the a^E+|ti=25)^B^n|w=4) 
and the c^E+|r;=35)~B^n|z;=12) vibrational levels, respectively (see text). 


3. Resonantly enhanced singlet-triplet mixing between B^II and 


3.1. Experimental procedure 

In order to further investigate the rotational, hne, and hyperhne structure of the two 
candidate state manifolds, we perform high resolution one-photon spectroscopy on 
^^Na^°K Feshbach molecules. Working with Feshbach molecules in contrast to free atoms 
enhances the coupling to the excited states, and deeply bound vibrational levels can be 
efficiently addressed. Instead of seconds, resonant excitation of Feshbach molecules 
typically occurs on a time scale of tens of microseconds. 

To create Feshbach molecules we use radiofrequency (rf) association as described 
in 20 . In short, after preparation of the mixture in the optical dipole trap, 
atoms are transferred to the |9/2, —7/2) state. Then, a uniform magnetic held 
of 85.7 G is applied to the atoms, close to a s-wave Feshbach resonance in the 
|/Na=l, /k=9/2, m/j^= —9/2) scattering channel caused by the closed channel 

molecular state a^S’''|n=19, iV=0, J=l, F=9/2, m^= —7/2). Here, F is the total 
molecular angular momentum including the Na and K nuclear spins, and mp its 
projection along the quantization axis. The quantization axis is set by the magnetic 
held in the vertical direction (z-axis). We associate Feshbach molecules by driving the 
free-bound transition from |/Na=l, /k=9/2, —7/2) to the molecular state 

causing the Feshbach resonance in |/Na = = 1) /k = 9/2,m/j^= —9/2), with the 

additional binding energy being provided by the rf photon. At the given magnetic held, 
we measure the binding energy of the Feshbach molecules to be 80 kHz and convert 











Normalized molecule number Q. Normalized molecule number c Normalized molecule number 


9 


J=1 J=Q J=2 J=2 

iZ_ ▼. T 



150 160 170 180 190 200 210 

Laser frequency (GHz) - 360 THz 


220 


230 


770 MHz 


550 MHz 



Figure 4. High resolution single-photon spectroscopy on Feshbach molecules near 
360 THz, exploring the coupled B^H |v=4)~c^E+ |ii=25)~b^no |i;=59) state manifold, 
(a) 90 GHz wide scan exposing a rich fine structure of nine loss features. The 
indicated loss features at 360.190 THz and 360.196 THz are dominantly c^E+ states 
with J=1 and J=2, respectively. Their hyperhne structure is shown in high resolutsion 
scans of panels (b) and (c). The two features at 360.221 THz and 360.229THz 
correspond to the dominantly singlet B^n|J=l) and B^n|J=2) level, respectively. 
Their Zeeman and hyperhne structure is shown in panels (d) and (e). In (b)-(e), the 
gray solid lines represent a calculation of the hyperhne structure, including the Zeeman 
interaction (see text). The blue, green and red lines indicate the spectroscopic weight 
of the individual contributing hyperhne features with mp= —9/2, —7/2 and —5/2, 
respectively, corresponding to tT“, tt and (t+ light. The gray shading in (d) marks 
the mj= -1-1 loss feature that is used for ground state spectroscopy (see Figs. and 
1^. The relevant sub-features (green) have quantum numbers mFi= — 1/2 (left) and 
mFi=l/2 (right). 
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Figure 5. High resolution single-photon spectroscopy near 372.5 THz, exploring 
the coupled B^H |ri=12)^c^E+|u=35)~b^no |n=66) state manifold, (a) A 100 GHz 
scan reveals eight loss features, (b) and (d) displays the hyperfine structure of the 
two J=1 states. Significant singlet-triplet mixing is evident due to the similarity of 
the observed features. The differing width of the hyperfine structure indicates the 
differing admixture of the hyperfine-rich c^S^ state. Likewise, the J=2 states shown 
in (c) and (e) share the same hyperfine features. Gray solid lines, and blue, green and 
red lines in (b)-(e) as in Fig.|^ The gray shading in (b) marks the sublevel that is used 
for dark state spectroscopy (see Fig. |^. The relevant sub-feature (green) has quantum 
numbers Fi=5/2 and mi?j=l/2. 
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about 10% of the atoms to yield about 10^ Feshbach molecules. 

The associated Feshbach molecules are illuminated by laser 1 (see Fig. to drive 
transitions to excited states of ^^Na^°K. Resonant coupling leads to a loss of Feshbach 
molecules. The intensity of the laser and the exposure time are chosen to clearly reveal 
the sub-structure of the excited states without saturating the transitions. The laser 
polarization is fixed to be diagonal with respect to the bias magnetic field. After 
exposure to the probe laser, an absorption image is taken and the number of remaining 
Feshbach molecules counted as a function of laser detuning. For absorption imaging, we 
use light resonant with the atomic transition. Absorption of a first photon breaks 
the weakly bound Feshbach molecules, followed by scattering of additional photons on 
the atomic transition. 

To achieve high spectral resolution required for precision spectroscopy, the 
frequency of the probe laser is locked to a high-finesse ultralow expansion (ULE) cavity 
in a master-slave configuration. The master laser is directly locked to a transmission 
mode of the ULE cavity by a Pound-Drever-Hall lock. The probe laser is phase-locked 
to the master laser with a variable frequency offset. For both lasers, we use home-built 
grating stabilized diode lasers. This configuration allows flexible adjustment of the probe 
laser detuning over a range of 10 GHz with kHz precision. The cavity has a finesse of 
~35,000 near 360.2 THz (832.2 nm), and ~15,000 near 372.6 THz (804.7 nm), suitable 
for addressing both candidate states. The temperature of the cavity is stabilized for 
long term frequency stability. Absolute frequency determination is performed for the 
most part with the help of a commercial wavemeter of 300 MHz accuracy as described 
above. One particular single-photon transition is determined to 3 MHz accuracy using 
a frequency comb, as discussed below. 


3.2. Analysis of the fine structure 

Figures 1^ and show the observed fine (a) and hyperfine (b)-(e) structure of the two 
state manifolds at 360 THz and 372 THz, respectively. In order to understand the 
spectral structure, we analyze the perturbations among nearly degenerate electronically 
excited states, following the approach of Ref. 53 . As a result of the theoretical analysis 


presented here, we are able to reproduce the observed fine and hyperfine structure, 
to characterize the electronic state decomposition, and to assign angular momentum 
quantum numbers. 

The fine structures of the 360 THz and 372 THz manifolds in Figs. |^a) and[^a) 
show nine and eight significant loss features, respectively. The state diagram of Fig. 
illustrates the various molecular states that are involved in the fine structure. Starting 
with non-rotating triplet Feshbach molecules in a A^=0, J=1 state of positive parity, 
only electronically excited states with negative parity and J=0, J=1 or J=2 can be 
reached. The c^S+ electronic state, being close to Hund’s case b), offers four such states: 
|A^=1, J=(0,1, 2)) and |A^=3, J=2), where the latter is coupled to |A^=1, J=2) via spin- 
spin coupling (vertical dotted line in Fig. [^. The B^H electronic state contributes two 
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Figure 6. Diagram of deperturbed states relevant for the chosen two-photon pathway 
connecting the Feshbach state in a^S+ to the singlet ground state X^E+. For the 
ground (electronically excited) state manifolds, only positive (negative) parity states 
are relevant and shown here. Accordingly, electronically excited states with odd (even) 
J have e (/) parity. B^II states are coupled via spin-orbit coupling (dashed lines) 
to c^E+, which in turn is coupled to b^IIo via spin-orbit coupling (short dashes), L- 
uncoupling (long dashes) or both (dot-dashed line). The vertical dotted line connecting 
c^E+ |fV=3, J=2) and |fV=l, J=2) represents spin-spin coupling. Transitions used for 
two-photon coupling to the J=0 absolute rovibrational ground state are indicated. 
To reach J=2 of the singlet ground state, one can either choose J=1 or J=2 in the 
intermediate excited state. 


states, the e state with J=1 and the / state with J=2. The remaining three (360 THz 
manifold) and two (372 THz manifold) features are associated with the b^Ho electronic 
state, introducing three states with J=0, J=l, and J=2. As only negative parity states 
are involved in the excited state, odd J corresponds to e parity and even J to f parity. 
The e// label will be omitted in the following discussion. 

The hne structure around 360 THz is associated with the perturbed complex 
B^n|n=4)~c^S+|r;=25)~b^no|n=59), while the 372THz manifold is caused by 
B^n|f=12)~c^S’''1^=35)~b^no|n=66). For clarity of this presentation, we neglect the 
influence of the other b^ni 2 states that are far off-resonant, a distance of 15cm“^ and 
30cm“^ higher in energy 55. However, they are included in the calculations that 


reproduce the observed spectra, where they cause minor shifts of lines. 

The Hamiltonian governing the hne structure of the perturbed B~c~ 

is 


b complex 


53,54,61 




The deperturbed energy levels F^B,c,b of electronic states and their rotational energies 
(with rotational constants i?B,c,b) are included in E-q and E^, are taken to 
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be the deperturbed energies in the absence of spin-orbit conpling of B^n|J=l) and 
c^E+|A^=l, J=l), respectively. i?b is the depertnrbed energy of the b^IIo |^= 0 ) state. 
Within we inclnde spin-spin conpling iifgg with conpling constant A, and the weak 

spin-rotation conpling with constant 7 , which is diagonal in the chosen Hnnd’s case 
b) basis. The conplings are schematically shown in Fig. States B^II and c^S’*' are 
connected via spin-orbit conpling (conpling constant .^bc)- Also the electronic states 
b^IIo and c^E+ are coupled via spin-orbit coupling H^q (coupling constant .^bc), and in 
addition, for J 7 ^ 0, via weak L-uncoupling (coupling constant Bl). The later has 
a minute influence for the low rotational states considered here, and will be neglected 
in the discussion, but is included in the calculation. 

All of the above perturbations conserve total angular momentum J, so that the full 
Hamiltonian can be broken down into sub-Hamiltonians that only act within sub-spaces 
with hxed J. In the basis {c^S’''|iV=l, J=0),b^Hol J=0)}, the Hamiltonian for J=0 is 

Ec - j -2\ -A/2^bc 

— V^^hc -Eb 

For J=l, neglecting L-uncoupling, another 2x2 matrix connects the states in the basis 
{B^n| J=l), c^S+|iV=l, J=l)} via spin-orbit coupling. 




Hj=,= 


Eb 'Cbc 

'bBc Ef, 


L-uncoupling causes weak coupling of c^S+|iV=l, J=l) to b^no|J=l), which is 
responsible for the extremely weak and narrow J=1 feature near 360.162 THz in 
Fig-i a). The corresponding feature in the 372 THz manifold was not found, presumably 
due to even weaker L-uncoupling. Finally, the J=2 sub-space is governed by a 
4x4 matrix in the basis {B^n| J=2), c^S+|iV=l, J=2), c^S’''|iV=3, J=2), b^Ho] J=2)}, 
neglecting influence of far off-resonant states b^Hi 2 , 


^ Eb + 4i?B 





y^^Bc 

Ec + 27 - |A 

2 v /6 \ 

\J\ihc 



Ec + lOBc - 37 

'6l 

be 


0 ^ 



Eh + QBh j 


As a starting point for the analysis of the spectra, we use deperturbed energy levels 


and rotational constants, obtained by mass-scaling of existing ^^Na^^K data on B^H 41 


c^S+ 52,57 and b^Ho (5^[5^[^ states. b^Holn = 66 ) has not been observed before and 
falls outside the validity of the Dunham expansion given in 57 . The above perturbations 


shift the energy levels, and the corresponding coupling constants serve as variables to 
match the theory with the experimental data. An initial estimate of ^bc = 0.58 cm“^ 
for the 360 THz manifold and 0.16 cm“^ for the 372 THz manifold was discussed in the 
previous section (57|. The value of .^bc = 0.77 for the 360 THz manifold is obtained from 


measured values in 57 after weighting with the proper Franck-Condon factors. Values 
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of 7 and Bl for can be fonnd in 52 and serve as initial values here. By varying 


the spin-spin coupling constant A of the unknown spin-orbit coupling ,^bc of the 

372 THz manifold, and the deperturbed energy levels £'B,c,b, excellent agreement with 
the observed hne and hyperhne structure is obtained (see Figs. |^and[^. 

The J quantum numbers of the observed hne structure features are identihed by 
considering their distinct hyperhne and Zeeman sub-structure. The two J=0 states 
are featureless, and J=1 and J=2 states show characteristic patterns, to be analyzed 
below. For either of the two manifolds, strong b~c spin-orbit coupling explains the 
unusual presence of a J=0 state in the midst of states with higher J. 

For the present task of hnding a two-photon pathway from Feshbach molecules to 
the absolute singlet ground state, the J=1 states are highly relevant. In the 360 THz 
manifold (see Fig. |^a)) the two J=1 features are well separated by 1 cm“^. This is 
signihcantly larger than the expected strength of spin-orbit coupling for these states 
given the initial estimate of ^bc = 0.16 cm“^ (our analysis yields 0.27cm“^, see below). 
Thus the states are relatively weakly mixed. As a consequence, the high resolution scans 
of the two J=1 features in Fig. |^b) and (d) show weak resemblance in their structure. 
Nevertheless, the very fact that two J=1 features are observable is a consequence 
of spin-orbit coupling, because a purely singlet state would not be detectable in our 
The analysis reveals that state (b) has 92% triplet character, and therefore 


scans. 


displays a characteristic triplet hyperhne structure; state (d) has 92% B^H character, 
and accordingly the most prominent feature is a Zeeman triplet of energy levels with 
magnetic quantum numbers mj= —1, 0, and -|-1. 

In stark contrast, the 372 THz manifold (see Fig. [^a)) contains two pairs of states, 
two J=1 states ((b) and (d)) and two J=2 states ((c) and (e)), whose respective 
hyperhne structure is nearly identical. This “line doubling” of an apparent c^E’*' state 
is well-known from spectroscopic studies of NaK 41,49-53 , and it results from strong 


mixing of c^E"*" with B^H. The energy diherence between the two J=1 features is about 
1.23 cm“^, comparable to twice the spin-orbit matrix element ^bc = 0.58. As the |J=1) 
states can approximately be considered to form a two-level system (see i7j=i above), 
we can write each of them as a superposition: 


I J=l)^^^ = a |c3e+, N=1, J=1) + f3 jB^H, J=l) (2) 

I J=l)^^^ = -13* |c3E+, N=1, J=1) + a* |B^n, J=l) (3) 

As the hyperhne structure is largely dominated by the c^E"*" contribution, the ratio 
of the widths of the hyperhne features reveals the value of |a|^/ \/3f ~ 1.6. The full 
perturbation analysis including hyperhne and Zeeman structure gives |a|^/ \(3f = 1.8, 
implying a = 0.80 and f3 = 0.60. We thus have identihed a pair of states with a nearly 
even, 64%-36% mixing between singlet and triplet states. 


3.3. Analysis of the hyperfine structure 

To understand the hyperhne and Zeeman structure found in the high-resolution scans 
of Figs. ib) -(e) and Figs. lli:b)-( e), we include in the total Hamiltonian the Zeeman 
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interaction of states B^II and and the hyperfine interaction within c^S+, 

consisting of the Fermi contact interaction of the electronic spin with the sodium and 


= -ff«„ + Hi + HI + + J?!;, K. 


the potassium nucleus (if^f and j^) as described in Ref. 


53 


The Zeeman Hamiltonian for the B^H state is 61 


Hi — ^bLzB, 


where hb is the Bohr magneton and the 2;-component of the molecule’s orbital angular 
momentum in the laboratory frame. is diagonal in the |Jmj) basis of rotational 
states of R^n|n) with matrix elements 

{Jmj\Hl\Jmj) = (4) 

describing the projection of the one unit of angular momentum along the internuclear 
axis onto the magnetic held in the laboratory 2;-axis; rrtj denotes the 2;-projection of J. 
In particular, the magnetic moment of a |J=l,mj=l) state is fiB/‘^, causing a Zeeman 
splitting of about 60 MHz, close to what we observe for the dominantly B^H state in 
Fig. I^d). Additional splitting and broadening is due to the admixture of c^S’*', which 
introduces hyperhne and additional Zeeman shifts. For the c^S+ state, the Zeeman 
term is 


= f^B9sS.B 


with Sz the ^-component of the electron spin the laboratory frame, and gs its ^f-factor. 
Note that is not diagonal in the Hund’s case b) basis, which mixes states of differing 
spin projection along the internuclear axis 61 1^. The hyperhne Hamiltonians are 

= ttNa^^Na ‘ S and H 


H\ 


hf,Na 


hf,K 


= oikIk ■ S, 


where /Na and Ik denote the nuclear spin for ^^Na and respectively. A major 
simplihcation of the analysis is ahorded by the large splitting between states of 
diherent J, compared to the hyperhne (and Zeeman) splitting (aNa ~ 340 MHz and 
ctK ~ 30 MHz). J therefore remains a good quantum number. 

The dominant contribution of the hyperhne interaction in c^S’*' originates from 
the ^^Na nucleus, which has a nuclear spin JNa=3/2. This leads to three (four) 
hyperhne manifolds for J=1 (J=2) at zero held, labelled by the quantum number Fi, 
corresponding to the addition of the sodium nuclear spin to the total angular momentum 
J excluding nuclear spins, Fi = J + I^g,. For J=l, possible values are Fi = 1/2, 3/2, 
and 5/2; for J=2 additionally Fi = 112 occurs. The appropriate hyperhne constants for 
various vibrational levels of c^E’*', including the ones relevant here (|n=25) and |n=35)) 
were measured in Refs. for ^^Na^^K. To include the hyperhne interaction of 

electron spins with the potassium nucleus, one introduces the total angular momentum 
F = Fi +Ik = J + Inii + Ik, with Ik = 4, causing new levels to appear from each Fi level 
(for example, Fi=l/2 leads to F=7/2 and F=9/2). Finally, the Zeeman interaction at 
85.7 G further splits these lines. The only remaining strictly “good” quantum number 
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Parameter 

B n=4)~c n=25)~b n=59) 

B n=12)~c n=35)~b n=66) 

Eb 

17289.158 

17701.427 

E, 

17288.338 

17701.074 

E^ 

17287.308 

17699.680 

Bb 

0.0662 

0.0551 

Bo 

0.0461 

0.0392 

5b 

0.0624 

0.0540 

Cbc 

0.27 

0.5899 (J=l), 0.5833 (J=2) 

^bc 

0.766 

0.68 

A 

-0.174 

-0.55 

7 

0.0022 

0.0215 

Bl 

0.0035 

0.006 

«Na 

0.0112 

0.0114 


-1.2 10-=^ 

-0.94 10-3 


Table 1. Molecular constants and interaction parameters for the two resonant 
B^n^c^S+^b^IIo manifolds. Calculated line positions of the fine structure agree with 
the experimental data to within the experimental accuracy of ^300 MHz, dominated 
by wavemeter error. All values are in units of cm“^. 


is the z-component of the total angular momentum, mp- Still, for well-separated levels, 
F or Fi and the ^-component mp^ are often approximately dehned. 

In order to understand the observed hyperhne spectra, we analyze which states are 
accessible starting from the initial Feshbach molecular state with mp= —7/2. Possible 
hnal states that can be reached via diagonal polarization are mp= —5/2, —"J/2, —9/2. 
An additional approximate selection rule comes from the fact that the closed-channel 
Feshbach molecular state is predominantly in the total angular momentum state F=9/2, 
the open-channel is in F=7/2, so that transitions to excited states F=13/2 or F=3/2 are 
forbidden. Working in the nuclear spin-decoupled basis of \N, J, mj, iNa, Ik, ’^/k) 
for c^S’*', and noticing that the Feshbach molecules have |iV=0, J=l), the line strength 
is described by the Honl-London factor {Jmj\lplm'j) p. The symbol in brackets 
denotes a Clebsch-Gordan matrix element, and p describes the laser polarization, where 
p = 0 stands for linear, and p = ±1 for circularly polarized light. The calculated 
spectra are shown above the measured data in Fig. (b)-(e) and Fig. (b)-(e). They 
reproduce the observed line structure well. Most lines consist of several individual 
hyperhne components, indicated by individual lines within the theoretical spectrum. 
Table summarizes the parameters used to reproduce the observed spectra. 

We can compare our results for the hyperhne constants with the literature on 
23Na39K. For c^S+|n=25) and |n=35). Ref. 52 found ONa = 312(6) MHz (in agreement 

For ^^Na^°K, our analysis yields ONa = 


with 


_|) and 342(6) MHz, respectively. 

334(4) MHz for |r;=25) and ONa = 342(2) MHz for |r;=35), the latter in excellent 
agreement with the work on ^^Na^^K 1^. The inhuence of the potassium nucleus on 
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the hyperfine structure had previously either not been observed 52 or only indirectly 


through broadening of Na-dominated hyperhne lines 53 , leading to estimates osqk < 
10 MHz and 18 MHz, respectively. From the known atomic hyperhne constant, the 


probability for hnding an s-electron near the nucleus was estimated to be 15% 53 
For this would yield = —22 MHz; note the negative sign due to the anomalous 
sign of the gf-factor of the nucleus. Here, we hnd a40K= — 37(4) MHz for |n=25) and 
a40K= — 28(4) MHz for |n=35). 

To summarize, we identihed promising intermediate states for two-photon coupling 
of Feshbach molecules to the absolute singlet ground state. In the next section we 
demonstrate this coupling directly. We hrst determine the previously uncertain location 
of the ground state via Autler-Townes spectroscopy and then perform coherent dark 
state spectroscopy for a precise measurement of the binding energy. 


4. Ground state spectroscopy of 

Before our study, the latest value on the position of the absolute rovibrational ground 
state of bosonic ^^Na^^K, relative to the dissociation limit of unbound atoms, was 
5211.75(10) cm“^ 44 . The relatively small error bar is a testament to the long history 


of high-resolution laser spectroscopy on NaK. For conventional heatpipe or molecular 
beam experiments, the main challenge is to precisely pin-point the dissociation energy, 
as the most weakly bound state of the triplet potential had not been directly observed. 
Working with Feshbach molecules, associated from free atoms, overcomes this problem. 
The measurement of interspecies Feshbach resonances between ^^Na and hxed the 
position of the weakest bound state relative to the dissociation energy 20,59 . Precise 


knowledge of the long-range van der Waals forces then allows in principle to assign 
absolute values to the binding energy of previously observed energy levels. 

Our observation of the singlet-rich feature in the B^n|u=4) c3S+|u=25)~b3n|u=59) 
complex (see Fig. |^d)) allows for a direct measurement of the singlet rovibrational 
ground state binding energy using two-photon spectroscopy. The dominant spin singlet 
character of this state results in strong coupling to the singlet ground state, which in 
turn induces a large Autler-Townes splitting of the excited state. 

For conducting two-photon spectroscopy, similar to the experimental procedure of 
the previous section, we hrst prepare an ultracold gas of ^^Na^°K Feshbach molecules 
in a crossed optical dipole trap. Then, the molecules are simultaneously exposed 
to a weak up-leg laser 1 resonantly addressing the singlet dominated J=1 state 
of the B^n|u=4)~c^S+|u=25) branch, and a strong down-leg laser 2 coupling the 
X^S'''|u=0, J=0) state to this intermediate state (see Fig. [^. After a set exposure 
time, an absorption image of the remaining Feshbach molecules is taken by imaging 
the component of the molecules. The number of remaining Feshbach molecules is 
recorded as a function of the down-leg laser detuning. 

In the absence of the down-leg light, the up-leg laser would simply remove most 
Feshbach molecules from the trap. However, when the down-leg laser is strong and 
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B^n ~c^i+ 



•-0 


B^n ~c^i+ 



Figure 7. Two-photon ground state spectroscopy of via the J=1 level at 

360.22115 THz (compare Fig. |^. (a) The up-leg frequency wi is kept on resonance, 
while the down-leg frequency UJ 2 is scanned, varying the detuning S 2 (see schematic). 
Autler-Townes splitting of the excited state protects the Feshbach molecules from being 
excited, when uj 2 approaches the ground state resonance. The full-width half maximum 
of the protection window is about 1.3 GHz. The solid line serves as a guide to the eye. 
(b) Autler-Townes splitting is mapped out by keeping the down-leg frequency 102 on 
resonance and scanning the detuning (5i of the up-leg frequency wi. The line shape is 
reproduced with a model based on the three-level master equation (solid line). The 
peak distance yields a Rabi coupling of H 2 = 27r x 123(5) MHz. 

sufficiently close to resonance, it shifts the excited state via the AC Stark effect by 
more than a linewidth, thereby “protecting” Feshbach molecules from decay. This gives 
the criterion for the width of the protection window: The AC Stark shift Vt \/252 due 
to the down-leg coupling laser at detuning 82 needs to exceed the linewidth T of the 
excited state, so the detuning should not exceed Vt^/T. If the down-leg Rabi frequency 
is signihcantly larger than the excited state linewidth, it opens up a wide protection 
window, facilitating the observation of the singlet ground state. In this so-called Autler- 
Townes (AT) regime, the protection of Feshbach molecules originates from two-level 
interaction between the intermediate state and the target rovibrational ground state, 
and requirements on phase coherence between the two Raman lasers are relaxed. 

As a starting point for the study, the binding energy of the rovibrational ground 
state of ^^Na^°K is obtained via mass-scaling the best value for ^^Na^^K. For the chosen 
intermediate state, the resonant down-leg frequency lies around 516 THz, corresponding 
to 580 nm, in the visible range. For generating this wavelength, we use a continuous 
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wave dye laser operating with Rhodamine 6G dye. The broad tunability of a dye laser 
together with a large output power of >1W has been instrumental in locating the 
rovibrational ground state. 

In Fig.g a) the normalized Feshbach molecule number after exposure to the Raman 
lasers is plotted as a function of down-leg laser detuning. As the laser frequency is 
tuned close to the energy difference between the intermediate state and the rovibrational 
ground state, we observe a large window of Antler-Townes protection of about 1.3 GHz. 
This signals the observation of the absolute ground state. The width of the protection 
window was on the order of the uncertainty in the ground state binding energy, 
facilitating the observation of the ground state. 

To accurately measure the Autler-Townes splitting and determine the down-leg 
Rabi coupling G 2 5 we £x the down-leg laser to be on resonance between the intermediate 
state and the rovibrational ground state and scan the up-leg laser frequency. In Fig. 
(b), the normalized number of remaining Feshbach molecules is plotted as a function 
of the up-leg detuning. From the measured splitting and known beam parameters, we 
obtain a normalized down-leg Rabi frequency of CI 2 = 27rx 40(10) kHz x 
indicating strong coupling. 

With the rovibrational ground state located, we accurately measure the binding 
energy of the X^E+|n=0, J=0) state via coherent dark state spectroscopy. In the limit 
of weak down-leg Rabi coupling, the transition probability to excite Feshbach molecules 
into the intermediate state is coherently canceled when the Raman lasers satisfy the two- 
photon resonance condition. In contrast to Autler-Townes protection, the protection of 
Feshbach molecules is now due to electromagnetically induced transparancy (FIT) |^ , 
which relies on the relative phase coherence between the two Raman lasers. 

The Raman laser coherence is maintained by locking both lasers 1 and 2 to a 
common ULE cavity using the Pound-Drever-Hall technique. The cavity is coated for a 
hnesse of ~35,000 in the vicinity of the up-leg and the down-leg transitions, and has a 
free spectral range of 1.5 GHz. For providing fast frequency feedback to the coupling dye 
laser, an electro-optic modulator is inserted into the laser cavity and driven by a high 
voltage, high bandwidth amplifier. Frequency tuning of the dye laser is achieved with an 
external acousto-optic modulator. As in the previous section, for addressing the up-leg 
transition, we use two grating stabilized diode lasers in a master-slave conhguration. 

In Fig.g a) and (b), we show the observed dark state resonance lineshapes, using the 
360 THz singlet dominanted J=1 feature in Fig. |^d), as well as the lower J=1 feature 
of the 372 THz manifold in Fig. |^b), respectively. For detection, we record the number 
of remaining Feshbach molecules after exposure to the Raman lasers as a function of 
up-leg laser detuning. The down-leg detuning is hxed to be on resonance. The observed 
lineshapes are htted by a simplihed three-level master equation model, which takes the 
down-leg Rabi frequency as a free parameter. For convenient absolute calibration of 
the down-leg frequency UJ 2 , we employ iodine spectroscopy and reference to the iodine 
atlas [^, cross-checked with the program Iodine Spec 66 . The frequency difference 


between the intermediate state feature of Fig. |^d) and the rovibrational ground state 
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a b 



Figure 8. Observation of the absolute ground state via two-photon dark state 
spectroscopy, (a) Dark state spectrum via the excited state feature at 360 THz for 
25 mW up-leg power (vertical polarization) and 20 fiW down-leg power (horizontal 
polarization). The dark state feature has a width of about 200 kHz. A three-level 
master equation model 64 matches the data for Rabi couplings Hi = 27rx 0.20(4) MHz, 
^2 = 27r X 2.6(5) MHz and a detuning ^2 = 27r x 2.5(3) MHz. This yields normalized 
Rabi couplings of Hi = 27r x 0.25(5) kHz x •^///(mW/cm^) and H 2 = Stt x 100(20) kHz 
X a///( mW/cm^). The linewidth of the excited state is Tie) = 27r x 10(1) MHz. (b) 
Dark state spectrum observed via the excited state feature at 372 THz using vertical 
(diagonal) polarization for the up-leg (down-leg) laser. The master-equation model 
yields the couplings Hi = 27r x 0.4(1) MHz, H 2 = 27r x 2.5(5) MHz and a detuning 
(52 = Stt X 0.8(1) MHz. Normalized Rabi couplings are Hi = 27r x 0.45(10) kHz 
X a///( mW/cm^) and H 2 = 27r x 25(5) kHz x a///( mW/cm^). The linewidth of the 
excited state is Tie) = 27r x 9(1) MHz. 


was measured to be a ;2 = 27r x 516.472613(3) THz, accurate to the MHz-level. The 
corresponding up-leg frequency tui = 27r x 360.221309(1) THz was calibrated using an 
optical frequency comb. We obtain a binding energy of Z1 q^^= 5212.0447(1) cm“^ for the 
rovibrational singlet ground state of ^^Na^°K, relative to the dissociation threshold and 
the hyperhne center of mass of the constitutent atoms. Adding the zero-point energy 
61.5710 cm“^, the depth of the ground state potential is Zle^^=5273.6157(1) cm“^. 
This corresponds to a thousand-fold reduction in the uncertainty compared to previous 


determinations 44 


5. Conclusion 

In this work, we have identihed a suitable pathway towards the production of ultracold 
ground state molecules of fermionic ^^Na^'^K, starting with Feshbach associated pairs 
of ultracold atoms. One-photon spectroscopy directly revealed the c^S+ triplet state of 
NaK, whose levels had only been observed before in the presence of strong perturbation 
by nearby singlet states. We identihed two strongly perturbed c^E+~B^n~b^no 
manifolds in ^^Na^°K, which feature levels of strongly mixed singlet and triplet character. 
Hyperhne-resolved spectroscopy revealed in one case almost equal singlet-triplet mixing. 
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These excited energy levels were shown to be ideal stepping stones for a two-photon 
process that connects the weakly bound Feshbach molecular state to the |n=0, J=0) 
absolute rovibrational singlet ground state. We performed Antler-Townes spectroscopy 
to detect the absolute ground state, and demonstrated coherent two-photon dark state 
spectroscopy for a precision measurement of the singlet ground state binding energy. 
The two resonantly mixed excited state manifolds should allow an efficient stimulated 
rapid adiabatic passage (STIRAP) from Feshbach molecules to absolute singlet ground 
state molecules. In this ground state, ^^Na^°K molecules are chemically stable against 
binary collisions 37 and possess a large induced electric dipole moment of 2.72 Debye. 
The results presented here pave the way towards the creation of a strongly interacting 


dipolar Fermi gas of chemically stable NaK molecules 67 
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